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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USEBS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 XJSEPS Group, in con- 
nection therewith. 
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1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 



Program No. Date 

Program Name: 



1. Does the abstract adequately describe what the program is and 

what it does? Yes No_ 

Comment 



2. Does the program do what the abstract says ? Yes No_ 

Comment 

3. Is the Description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient 

detail? Yes No_ 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently 

understandable ? Yes No_ 

Comment 

5. Does the source program compile satisfactorily (if applicable) ? Yes No_ 

Comment 

6. Does the object program run satisfactorily? Yes No_ 

Comment 

7. Number of test cases run 



Mr. Robert J. Robinson (PREP) Company 
Marquette University 

Computing Center Address 
1515 W. Wisconsin Avenue 

Milwaukee 3, Wisconsin User Group 

Code 



Are any restrictions as to data, size, range, etc. covered 

adequately in description? Yes No_ 

Comment 

8. Does the Program meet the minimal standards of the 1620 Users 

Group? Yes No_ 

Comment 

9. Please list any suggestions to improve the usefulness of the 
program. These will be passed on to the author for his consideration. 
Comment 

Please return to: Your Name 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
TO PARTICIPATE IN THIS EVALUATION. 



♦ ♦ 
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Kansas State Teachers College 
Emporia, Kansas 

Title ; Linear Programming 

Subject Classifica ti on: 10.1. 013 

Author ; Org anization : Roy Gallup; Kansas State Teachers College 

Di rect Inquiries t o: Roy Gallup 

KSTC Computer Center 
Emporia, Kansas, 66301 

Phone: DI2-5000 EX. 222 

Descriptio n: This program solves Linear Programming, problems 

using the simplex algorithm. All the identities, 
artificial variables and necessary changes to 
obtain a optimal feasible solution are added by 
the program. All that the user has to do is enter 
the coefficients of the variables a^ ^ , the cost 
coefficients c^, and the requirements bj_ . The 
optimal feasible solution is punched out along 
with the cost, basis vectors and the basic variable 

Restrictions : The maximum size problem that can be solved is; 

on a 20K machine, 13 equations and 13 variables, 
on a 40K machine, 31 equations and 31 variables , and 
on a 60K machine, 40 equations and 40 variables. 



Equipment Specifications : 



1 . Card system. 

2. Automatic divide. 

3. Indirect addressing. 

4. Floating point hardware . 

5. TNF and TNS instructions. 

6 . Language ; SPS . 

7. Subroutines; none. 

Running Time : A problem with 4 equations and 4 variables 

takes' approximately 5 seconds per iteration. 
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I. PROGRAM DESCRIPTION. 

This program solves linear programming problems using the 
simplex algorithm in which all identities and artificial 
variables are added by the program. The only thing the 
user has to do is to enter the coefficients of the variables 
aj_j, the cost coefficients Cj, and the requirements b^. 

The optimal feasible solution is punched out along with the 
cost, basis vectors and basic variables. 

The mathematical formulation of a linear programming 
problem is: 

Find values for X-j, X p , ... X which satisfy the 
following conditions; n 

a 11 X 1 + a 12 X 2 + ... a 1n X n ^ b 1 

a 21 X 1 * a 22 X 2 + a 2n X n % b 2 

a m1 X 1 + a m2 X 2 + ••• a mrPn # b m 

where Xj for all values of j and 

b^ > for all values of i. 

max z = c-jX^ + C 2 X 2 + c n X n 

where a^ , b i , and are constants. 

A. Machine Requirements. 

1 . Card input and output. 

-2. Floating point hardware. 

3. Automatic divide hardware. 

v 4. Indirect addressing. 

-5. TWF and TNS instructions. 

-6. Source language: SPS. 

("■ 7 ♦ Subroutines: none. 

8 . All computations are in floating point using a 
8 digit mantissa. 
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' 9'. Any size storage can be used. The size of the problem 

is restricted by the following relationship: 

(2 y j + i + l)(j + 3) memory - 7300 

r - 10 
where: i is the number of variables. 

j is the number of equations including 

the cost row. 
memory is 20,000, 40,000, or 60,000. 

B. Error Stops and Procedures. 

1. If the parameter card is missing, the message "FIRST 
CARD IS MISSING-" will be typed out and the computer 
will halt. To restart — lift the remaining cards out 
of the read hopper, depress non-process run out to run 
cards out of reader. Insert parameter card in front 
of the data and place it in read hopper. Push START 
on the console and push READER START on 1622. 

2. If a card is missing within a equation or the requirement 
card is out of order or missing the message "INCORRECT 
NUMBER OF EQUATIONS" will be typed out and the computer 
will halt. To restart — lift the remaining cards out 

of the read hopper, depress non-process run out to run 
cards out of the reader. Correct the data cards and 
place them in the read hopper. Push START on the 
console and push READER START on 1622. 

C. To Change the Size of the Matrix. 

To change the size. of the matrix when using a 1620 
with more than 20,000 storage positions, columns 10-11 
in the card number 3 (in column 80), which contains 
j (the number of equations including the cost row), 
should have 18 for 20,000; 31 for 40,000 and 40 for 
60,000 storage positions. It should be noted that 
the size includes space for the identity and artificial 
variables . 

(' - D . Execution Time. 

The time depends on the size of the problem to be 
solved. A problem with 3 variables and 4 equations 
and two artificial variables takes approximately 
3 to 5 seconds per iteration. 
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II. THE PROGRAM. 

A. Loading and Operating Instructions. 

1. Set PROGRAM switches to STOP. 

2. Set console switch 1; 

ON: Punch out complete matrix for each tableau or 

iteration. 
OFF : Supress punching . 

3. Set console switch 2; 

ON: Type value of max z for each tableau or iteration. 
OFF: Supress typing. 

4. Place blank cards in punch hopper. 

5. Place object deck and data cards in read hopper. 

6. Depress LOAD button on 1622 (program has a self- 
clearing routine) 

7. After object deck is loaded, the computer will halt. 
To begin execution of the program push START on the 
console . 

8. After all the data has been read in and the computations 
completed, the answers will be punched out. If the 
solution is unbounded or degenerate, a message card is 
punched out with "NO FEASIBLE SOLUTION" on it. The 
program is ready for another set of data. 



B. Data Preparation. 

1 . The first card is the parameter card. 

col. 1-2 contains the number of variables (i). 

col. 3-4 contains the number of equations including 

the cost row. (J), 
col. 50-80 can be used for identification purposes. 

2. The remainder of the data cards contain (one per card) 
the coefficients of the variables, the coefficients 

of the cost variables and the coefficients of the 
requirements in E14.8 format, i.e., +xxxxxxxxE+yy 
where the x's are the mantissa and the y's are the 
exponent. 
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The card format for the coefficients of the variables 
and the coefficients of the costs variables is; 

col. 1 contains the sign of the mantissa, 

col. 2 contains a decimal point, 

col. 3-10 contains the mantissa, 

col. 11 contains the E. 

12 contains the sign of the exponent. 



col . 
col . 



13-14 contains the exponent. 



Examples ; 

+62.5 

-62.5 
+ .732 

+ .073 
-.0053 



would be 
would be 
would be 
would be 
would be 



+.62500000E+02 
- . 62500000E+02 
+.73200000E+00 
+.73000000E-01 
-.53000000E-02 



The card format for the requirements is the same 
with one exception. The equality or inequality 
sign is placed in column one (1) and the rest of 
the data is shifted one place to the right. The 
equality consist of the equal sign "=" , which 
is denoted by a H =" sign or a 3-8 punch. The 
inequality consist of the less than equal sign 
" -=z ,f , which is denoted by a open parenthesis M (", 
or a 0-4-8 punch and the greater than equal sign 
" — which is denoted by a closed parenthesis 
")", or a 12-4-8 punch. 

col. 1 contains the equality sign. 

col. 2 contains the sign of the mantissa. 

col. 3 contains the decimal point. 

col. 4-11 contains the mantissa. 

col. 12 contains the E. 

col. 13 contains the sign of the exponent, 

col. 14-15 contains the exponent. 



Examples ; 



=+60.5 would be =+ . 60500000E+02 

2T-60.5 would be ) - . 60500000E+02 

■£+.065 would be ( + . 65000000E-01 
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3. Sample Equations. 



3X 1 + 5Xg £ 15 
-5X 1 - 2X g > -10* 

max z = 5X 1 + 3X ? 



5X + 10X ^ 50 



x 1 + 



x 2 ^ 1 
X 2 ^4 



min z = -X^ - X, 



u 

TT 



0203 

+.30000000E+01 

+.50000000E+01 

(+. 15000000E+02 

+.50000000E+01 

+.20000000E+01 

(+.10000000E+02 

+.50000000E+01 

+.30000000E+01 



0204 

+.50000000E+01 
+ . 10000000E+02 
(+.50000000E+02 
+.10000000E+01 
+. 10000000E+01 
)+.10000000E+01 
+.00000000E+01 
+. 10000000E+01 
(+.40000000E+01 
+. 10000000E+01 
+ . 10000000E+01 



Note all requirements must be positive, i.e. 



-5X, - 2X ^ -10 



would be 



5X 1 + 2X 2 ^: 10 



1 —2 

^Note that min z = max(-z) and therefore minz 
has to be multiplied thru by a -1 to give 



-X. 



max z = X 1 + x 



k 2* 



4. The total number of data cards for a problem can 
be determined by the following equation. 

parameter card + (i)(j) + j - 1 

where i is the number of variables. 

j is the number of equations including 
the cost row. 



- X 



Interpretation of the results. 

1 . When console switch 2 is on the MAX Z for each tableau 
or iteration is typed out except the last one which is 
punched out as the final answer. For example; 

TABLEAU 01 MAX Z = +. 00000000E+01 
TABLEAU 02 MAX Z = +. 1 0000000E+02 



is typed out while the third tableau is punched out. 
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2. Interpretation of a tableau when punch out in the 
complete matrix is; 



TABLEAU 01 



NUMBER GF TABLEAUX OR ITERATIONS 



+.00OOOOOOE+O1 
+ . ooooooonr+o 1 
+ .nooononnr+n l 
+•^0000000^+0 1 
+.30000000F+0 1 
+.OO00O000E+O1 
+.00000000E+01 
ROW 02 • • • • • 



• • • 



NUMBER OF THE ROW IN THE TABLEAU 



COST ROW COFFF I C I FNTS • 



+.00000O00E+O1 
+. onnOOOOOF+0 1 

+. noooooo^F+o i 
+ . i oooooooF+n i 



• • • 



BASIS VECTORS* 



+.P0000000E+01 • 
+.30000000E+01 • 

+ .40000000E+0 1 ••• 

ROW 03 

+ .OOOO000OF + 01 

+.30000000^+01 • 
+. 1 ^n0000OF+o? , 

+ . 30000000F+0 1 ... COEFFICIENTS OF FIRST FQUAT I ON. 

+.^ooonoooF+o l . 

+. 1 000OOOOE+O 1 • 

+ .00000000E + 01 • 

ROW 04 

+. OOOOOOOOE+O 1 

+. 40000000F+O 1 , 

+. 1 onnoonoF+n? • COFFF I C I FNTS OF FIRST FOUATION. 

+. ^OOOOOOOF+O 1 PIVOT POINT ... INDICATES PIVOT POINT FOR NEXT 
+. POOOOOOOF+0 1 • ITERATION. 

+.00000000E+01 • 

+ . 10000000E + 01 

ROW 05 

+. OOOOOOOOE+O l..... 

+.00000000E+01 • 
+. ooonooooF+o l • 

-.^OOOOOOOF+0 1 ... LAST ROW IS 7AJ) - C(J) FOR EAC 

30000000F+0 1 • COLUMN. 
+. OOOOOOOOE+O 1 
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3. The final answers that are punched out are interpreted 
as follows. 

TABLEAU 03 Number of tableaux 

or iterations 

COST BASIS VECTORS BASIC VARIABLES — heading 

.30000000E+01 2 +. 2368421 OE+01 

•50000000E+01 1 +.1 052631 6E+01 

MAX Z = +.1 236842 1E+02 

The summation of the COST times the BASIC VARIABLES 
equals MAX Z. 



4. Example of how to put matrix punch out into tableau 
form. Elements are taken from preceding page. 



TABLEAU 1 



ROW 


1 








5 


3 








ROW 


2 


CB 


VB 


b 1 


x 1 


x 2 


X 3 


X 4 


ROW 


3 





3 


15 


3 


5 


1 





ROW 


4 


o 


4 


10 


df 


2 





1 


ROW 


5 









-5 











'Circled element indicates pivot point. 
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p?03 ^AMPLF INPUT 1 

+. 30000000F+0 1 

+.50000000E+01 

(+. 15000000E+02 

+.50000000E+01 

+.20000000E+01 

< + . 10000000E+02 

+.50000000E+01 

+.3oooocoof+oi 



TABLEAU 01 
+.00000000E+01 

+.oooooooof+01 
+ .oor>oooonF+oi 

+ .F50000000F+01 
+.3000000CE+01 
+.00000000E+01 
+.00OO0O00E+01 

ROW 02 

+.00000000E+01 
+.00OO00OOE+01 

+.oooooooof+oi 
+. 1 ooooooof+o 1 

+.20000000E+01 
+.30000000E+01 
+.40000000E+0 1 
ROW 03 

+.0C000000E+01 
+ . 3<~>00000OE+0 1 
+ . 1 loooooof+o? 
+.3000PCOOF+01 
+ .P50000000E+0 1 
+. 10000000E+01 
+.OO0OO00OE+O1 
ROW 04 

+.00OOO00OE+O1 
+.40000000F+P1 
+ , 1 nnOOOOOF+02 
+ . B 50 0E+01 
+.20000P00F+01 
+.0OOOOOO0E+0 1 
+. lOOOOOOOE+Oi 



SAMPLE OUTPUT 1 
(MATRIX PUNCHOUT) 



PIVOT POINT 
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ROW 05 

+.00OOOO00E+O1 

+.00000000E+01 

+.00000000E+01 

-.^OOOOOOOF+O 1 

-.30000000F+01 

+.000000O0E+01 

+.OOOOOOOOF+01 

TABLEAU 02 

+.00OOO000E+01 

+.00000000E+0 1 

+.O0000000E+01 

+. LOOOOOOOE+OI 

+.3 n 000000E+01 

+.00000000E+01 

+.OO000O00E+O1 

ROW 02 

+.OO0OOOO0E+01 
+.0000OO0OE+0 1 
+.0OOO0000E+01 
+, 1 OOOOOOOE+C 1 
+.20000000E+01 
+.3O000000E+01 
+.40000000E+01 
ROW 03 

+.0OO0O0O0E+0 1 
+.30000000E+0 1 
+.90000000E+01 
+.OO00OO00E+01 

+.38000000E+01 PIVOT POINT 
+. 10000000E+01 
-.60000000E+00 
ROW 04 

+.50000000E+01 
+ . 10000000E+01 
+.20000000E+0 1 
+. 10000000E+01 
+. 4O000OO0F+00 
+.0OOO000OE+01 
+.20000000E+00 
ROW 05 

+.00000000E+0 1 
+.00000000E+0 1 
+ . 10000000E+02 
+.O000O0O0E+01 
-. 10000000E+01 
+.00000000E+01 
+. 10000000E+01 
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TABLEAU 03 
+.00000000E+01 
+. 00000000E+0 1 
+.OO0000O0E+O 1 
+.50000000E+01 
+.30000000E+01 
+.OOOOO00OE+01 
+.00000000E+D 1 
ROW 02 

+.OOOOOOOOE+0 1 
+.00000000E+01 
+.OOOOO0OOE+01 
■»•• 1O000O0GE+01 
+.20000000E+01 
+ . 3000000CF+0 1 
+.40000000F+01 
ROW 03 

+.30000000E+01 
+.20000000E+01 
+.23684210E+01 
+.00000000E+01 
+. 1OO00000E+01 
+.26315789E+00 
-. 15789473E-00 
ROW 04 

+.50000000E+01 
+ . 1000000CE+01 
+. 1 05263 16E+01 
+* 10000000E+01 
+.00000000E+01 
-.10526315E+00 
+.26315789E+00 
ROW 05 

+.O0OOO0O0E+01 
+.00000000E+01 
+ . 12368421E+02 
+.00000000E+0 1 
+.00000000E+01 
+.26315792E+00 
+.84210530E+00 



PRICE ' BASIS 
+.30000000E+O1 
+.50000000E+0 1 

MAX Z 



SAMPLE OUTPUT 1 (ANSWERS) 

TABLEAU 03 
VECTORS BASIC VARIABLES 

2 +.23684210E+01 
1 +. 10526316F+01 

: +. 12368421E+02 
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0204 SAMPLE INPUT 

+.50000000E+01 

+. 10000000E+02 

(+.50000000E+02 

+. 10000000E+01 

+. 1O0O000OE+01 

)+. iooooonOE+ei 

+.nnooonooE+oi 

+. 10000000E+01 
(+.40000000E+01 
+.1 OOOOOOOE+O 1 
+. 10000000E+01 



TABLEAU 01 

+.00000000E+01 

+.00000000E+01 

+.00000000E+01 

+. 1 OOOOOOOE+O 1 

+. 1 OOOOOOOE+O 1 

+.00000000E+01 

+.00000000E+01 

+.00OOOO0OE+01 

-. 10000000E+50 

ROW 02 

+.00O0O000E+O1 
+.00000000E+01 
+.000OOOOPE+O1 
+. 1OO00000E+01 
+.20000000E+01 
+.30000000E+01 
+.40000000E+01 
+. 50000000E+0 1 
+.60000000E+01 
ROW 03 

+ .00000000E+0 1 
+.30000000E+01 
+.50000000E+02 
+.50000000E+0 1 
+ . 10000000E+02 
+. 10000000E+01 
+.00000000E+01 
+.00000000E+01 
+.0OO0OOOOE+O1 



SAMPLE OUTPUT 2 
(MATRIX PUNCHOUT) 



o 
U 



u C 

H .2 

(U o 



o 
U 



2 « 

^ W z 



2 



????????? ???$?????■ ????°?°?? 

ggggg gg§§§I§gg ggggggggg gg|gggggg S § 

2§ 






sgggggggggsgggggggggisgggggggggs 

OClfiOOOOOC CCOOOOOCC OOOCCCCOO OCCCCOOOO_JCCCOOOOOC 
30fHtO(l-l(IOI('S----HO"C-SOI()<tOHCC«OJOC«CCO-0-IOCC--000» 
O • ••••••••0*********0*«*******0*********< ••••••••• 

a + + + + + + ++ i a + + + + + + i + + a + + ++ + + + + + a + + + + + +i + + 1- + +.+ + + + + + 1 




2 8 



f.ccoooocc ccooooccc ooSSSc?oo occocoooS CCOOOOOCC 

gggggr 





SgggggggggcgggggggggggggggggggSggggggggg^ggggggggg 

CCCOCOOCC COOCOOCCC OCOOCCCCO_JOCCCCCOOO CCOOOOCCC 

3 - c- ^- o-c-som^o-oc-csoo - ~ ~ c ~ o o m oc c - -ooo-jcco-«in«ii;<i 
o •••••••••o««««*««««o*********<*********o ********* 

ai+ + + + +i++a + + + + + + + + +a + +i 1 i+ + + + k+ + + + + + + +iqt + + + + + + + + + 
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ROW 02 

+.OO0000O0E+0 1 
+.00000000E+01 
+.00000000E+0 1 
+• 10000000E+01 
+.20000000E+0 1 
+.30000000E+0 1 
+.40000000E+01 
+.50000000E+0 1 
+.60000000E+0 1 
ROW 03 

+.OOOOOOOOE+01 
+.40000000E+0 1 
+.90000000E+0 1 
+.00000000E+01 
+ . 1 OOOOOOOE+O 1 
+.20000000E+00 
+. 10000000E+01 
+.00000000E+0 1 
1 OOOOOOOE+O 1 
ROW 04 

+• 1 OOOOOOOE+O 1 
+. 1 OOOOOOOE+O 1 
+• 10000000E+02 
+• 10000000E+0 1 
+.20000000E+0 1 
+.20000000E+00 
+.00000000E+01 
+.0O0OOO00E+0 1 
+.00000000E+0 1 
ROW 05 

+.00000000E+01 
+.50000000E+01 
+.40000000E+0 1 
+.00000000E+0 1 
+• 10000000E+01 
+.00000000E+0 1 
+.00000000E+01 
+♦ 10000000E+01 
+.0O000O0OE+01 
ROW 06 

+.OO0OOOO0E+O 1 
+.000O0OO0E+01 
+• 10000000E+02 
+.00000000E+0 1 
+• 10000000E+0 1 
+.20000000E+00 
+.00000000E+01 
+.0OOOOC0OE+01 
+• 10000000E+50 
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SAMPLE OUTPUT 2 (ANSWERS) 



TABLEAU 03 

PRICE BASIS VECTORS BASIC VARIABLES 

+.OOOOOCO0E+O1 4 +.90000000E+01 

+.10000000E+0I 1 +. 10000000E+02 

+.00OOOOO0E+0! 5 + . 40000000E+0 1 

MAX Z = +. 10000000E+02 



© 
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c 


c 1 


c 2 




C J 




°n 


-M 




-M 


C B 


VB 


b 


X 1 


X 2 




X J 




X 

n 


X . . 
n+1 




X 
n+s 


C B1 


X 1 


b 10 


Y 11 


Y 12 




Y U 




Y m 


Y 1 ,n+1 




Y 1 ,n+s 


C B2 


x 2 


b 20 


Y 21 


Y 22 




Y 2J 




Y 2n 


Y 2,n+1 




Y 

2, n+s 




























X 

m 


b mO 


Y m1 


Y m2 




Y mj 




Y 
mn 






Y 

m,n+s 






Z 


z r c f 


z 2~°2 




z r c j 




z n" c n 









where; j is the number of variables. 

n is the number of identities. 

s is the number of artificial .variables . 

mis the number of equations . 

The first column of the tableau (gB), gives the costs which 
correspond to the vectors in the column. 

The second column (VB), tells which vectors are In the basis. 
The third column gives the current value of th© requirement 
together with the value of the objective function for the basic 
feasible solution described by the given tableau. 
The remainder of the values are the Y for all vectors in X 
including any artificial variables which may have been adde9 +s 

The first row of the tableau gives the costs associated with the 
vectors. The last entry in each of these columns gives Zj-04 
for each vector. J 
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Q address 
determines 

the size of 
the matrix 



ADDSET 



seta up 

variable 

address 



SADDj 



sets up 
address 
3f work 



ar* a 



read in par- 
ameter card 
i and .1 



compute number 
of data cards 
to be read in 
(i)(.i)+( t i)-l 



sets up 
various 
counters 



set matrilx 
to zero's 



i - col\imns 1-2 Number of variables. 

j - columns 3-U Number of equations 
including cost 
equations . 



CC W 
set up number 
of variables 
in matrix 

x l> x 2> ••• X n-».s 








■MPPpiiPfPI 
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read in data 
cards 



NO 



store variable 
coefficient in 
right place in 
tableau 



does 

element belol 
in basis, 
^quireinent) 



YES 



B 



set up 
identity- 
element 




set identity 
element 
to minus one 



set flag 
at 
SETTAB 



strip out 
mantissa 
& exponent 




clear 
flag at 
SETTAB 




AA V 



store requlr- 
ment in ri^ht 
place in tableau 



NO 
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store then 
in first 
row of mat rix 



set up address 
to test 
identity- 
matrix 



YES 




set up ar 
variable and 
positive i 
matrix 



tif ica al 
and 

identity 



2 



set up 
column 
CB & VB 




12 H> 



P ! 

set up c 
to corpi 
row 


iddress 
ite 


MODI \ 


/ 


compute each 
element in 
row 

,1 J 
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ON — OUT 
switcK^v J L f 5 
1 P 



DPCj 




set a( 
to del 
pivot 


idress 
bermine 
column 


TEST \ 


1 


find smallest 
value in row 
z. - c. 


A10 \ 


i 




YES 



add 1 to 
negative 
count 
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7 




punch 
NO FEASIBIE- 
SOIUTION 




determine the 
smallest ele- 
ment in the 
work area 



SW1 




OUT 



SW22 



set address 
to type ou|t \>' 
max z 



set flag 
at 
SW2 




OUT 



set address 
to compute 
new tableau 



t 



divide 
pivot row 
by pivot point 



compute 
new 

tableau 
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OUT 
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edit floating 
point output 
to ElM.8 
format 




move max j 
to type- 
writer outnut 



write 
out 
max z 



move post 
to output 
area 



<3 



move max 2 
to card 



output area 




move basic 
Ojvariable tk 
output are 1 



set flag 
at 

PRICE 



punch 
out 
element 



NO 




increment 
row count 



' 4 > 





clear 1 




flag 4 




D0NS1 | 



)rancm 
back t< 
h02 



OUT 
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>uhch out num- \ 
ber of tableaux J 
; or iterajions J 



/punch oul 
(heading 



' "T 

, \ 


: — 7—^ 

1 


set add 
to edit 
variabl 


ress 
basic 

e 




XX 




clear fla^ 
at PRICE 


> 


move basis 
vector to 
output are 









punch 
output 
card 




set flag 
D0NE1 



OUT 

HO 
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LINEAR PROGRAM SIMPLEX ALGORITHM 
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ADD , ADD4 
CTPAR, I 
CTPAR,J1 
CTPAR, VARCT 
ADD7,T940 
AAA+1 1 ,TAB+20 
CT1 ,0. 10 
ADD5.T920 
WORK 1+9, FZERO 
ADD 6, ADD 7 

WORK+9, ADD5« 1 1 , COMPUTES ROW Z(J)-C(J> 
FMUL WORK+9 , ADD6 • 1 1 
FADD WORK 1+9, WORK+9 
A ADD5.S460 
A ADD6,S460 
CT1 , 1 , 10 
CT1 , Jl 
MOD 1 

WORK 1 +9 , TAB+20 , 7 
ADD4 , WORK 1+9,6 
AAA+1 1 , 10 
ADD7, 10 
ADD4,10 
CT2 , 1 , 10 
CT2, CTPAR 
AGAIN 
CT1 , 1 , 10 
ADD, 8 
XX 

*+24,ADD, 1 1 , TEST 
XX.,, SEE IF 

ADD , 1 
CT1 , 1 ♦ 10 

CT1 , CTPAR 
XX 

*+24,XX 

DPC 

PPADD+1 1 , 
SW1+12 
DONE, XX 

CT1 ,0, 10, DETERMINES PIVOT COLUMN 
CTPAR, 2, 10, ADD CONTAINS ADDRESS OF 
J , S460 , , SMALLEST VALUE IN ZJ-CJ. 
95 

99.T490 
ADD , 99 
ADD 1 ,99 
ADD 1,10 

WORK+9, ADD1 , 1 1 
FSUB WORK+9, ADD, 1 1 
BNF *+24,W0RK+7 



TF 
TF 
A 
A 

TF 

TFM 

TFM 

TF 

TFL 

TF 

TFL 



AM 

C 

BNE 

FSUB 

TFL 

AM 

AM 

AM 

AM 

C 

BNE 

TFM 

AM 

CF 

BNF 

SF 

AM 

AM 

C 

BNE 
BNF 

B 

TFM 

BC1 

BNF 

TFM 

SM 

M 

SF 
A 

TF 
TF 
AM 
TFL 



ROW ZJ-CJ TO 
ITS ALL POSITIVE 
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0?9?2 
02934 
02946 
02958 
02970 
02982 
02994 
03006 
03018 
03030 
03042 
03054 
03066 
03078 
03090 
03102 
031 14 
03126 
03138 
03150 
03162 
03 1 74 
031 86 
03198 
03210 
03222 
03234 
03246 
0325P 
03270 
03282 
03294 
03306 
03318 
03330 
03342 
03354 
03366 
03378 
03390 
0340? 
03414 
03426 
03438 
03450 
03462 
03474 
03486 
0349R 
03510 
03522 
03534 
03546 



26 06488 

I 1 06478 
24 06478 
47 02874 
16 064O? 
16 06478 
16 06494 

23 06476 
3? OO0Q5 
22 06488 
26 07127 
22 07127 
26 071-57 
26 07152 
26 07142 
33 03342 
26 07152 
33 03138 
06 071 17 
06 07107 

43 03186 
32 03138 

44 03222 

I I 06492 
32 03138 

43 03246 
32 03138 

44 03282 
1 1 06494 
32 03136 
44 03318 
06 0714K 
49 03354 
09 071 1 7 

32 03342 
06 0714K 
11 07142 
21 07152 
21 06488 

I 1 06478 
?4 0647R 
47 03126 
44 06124 

24 06492 
46 06124 
24 06494 
46 06124 

33 03534 
16 06478 
26 07122 
26 07147 

II 07147 
06 071 1 7 



07122 
000-1 
06464 
01200 
000-0 
000-0 
000-0 
07157 
Ooooo 

00099 
06488 
07157 
064 88 
07217 
00893 
00000 
07?17 
00000 
0715K 
0648Q 
071 08 
00000 
071 15 
000-1 
00000 
07098 
00000 
07105 
000-1 
00000 
03138 
06421 
00000 
071 07 
00000 
071 17 
-0010 
07157 
07157 
000-1 
06476 
01200 
03342 
06476 
01200 
06476 
01200 
000 00 
000-0 
00893 
07122 
-0010 
071 4P 



A 1 1 
0*1 1 " 



TF 
AM 

C 

BNE 
TFM 
TFM 
TFM 
M 

SF 

Q 

TF 

S 

TF 

TF 

TF 

CF 

TF 

CF 

TFL 

TFL 

BD 

SF 

8NF 

AM 

SF 

BD 

SF 

BNF 

AM 

SF 

BNF 

TFL 

R 

FDI V 

SF 

TFL 

AM 

A 

A 

AM 
C 

BNE 
BNF 

C 

BE 
C 

BE 

CF 

TFM 

TF 

TF 

AM 

TFL 



ADD. ADD 1 
CT1 . 1 , 10 
CT1 ,CTPAR 

TEST 

NEGCT.O. 10. DETERMINES PIVOT POINT 
CT1 .0. 10.ADP1 CONTAINS ADDRESS 
NEGC Tl.0.10, OF PIVOT POINT 
Jl ,S460 

95 

ADD . 99 
ADD2 t ADD 

ADD2.S460. .CONTAINS AODQESS OF a 

ADD4 . ADD 

ADD7.T940 

ADD5.SADD+1 l 

Al 1 

ADD7.T940 
A10 

WORK+9 . ADD7 . 1 1 
WORK 1 +9, ADD . 1 1 

* + 24 , WORK 
A10 

*+36.W0RK+7 
NEGC T . 1.10 
A10 

* + 24 .W0RK1 
A10 

*+36. WORK 1+7 
NFGCT1 .1.10 
A 1 

*+36, A10 

ADD 5 . FZERO . 6 

OMIT 

W0RK+O,W0RKl+9 
Al 1 

ADD5, WORK+9. 6 

ADDS, 1 

ADD7 , S460 

ADD.S460 

CT1 , 1 . 10 

CT1 , Jl 

A10-12 

ERR . A 1 1 

NEGC T . J 1 

ERR 

NEGC T 1 , Jl 

ERR 

A15 

CT1 .0. 10 
ADD 1 .SADD+l 1 
ADD6 . ADD 1 
ADD6. 1 
WORK+9, AD06, 1 1 
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03558 
03570 
03582 
03594 
03606 
03618 
03630 
03642 
03654 
03666 
03678 
03690 
03702 
03714 
03726 
03738 
03750 
03762 
03774 
03786 
03798 
0381 
03822 
03834 

C3858 
03870 
03882 
03RO4 
03906 
03918 
03Q30 
03Q4? 
03954 
03966 
03978 

040C2 
04C14 
04026 
04038 
04050 
04062 
04074 
04086 
04 098 
041 10 
04122 
041 34 
04146 
04158 
04 170 
04 182 



14 0714P 

46 03678 
14 0712K 

47 03630 
26 07122 
49 03678 
02 071 17 
32 03534 
44 03678 
26 07122 
1 1 06478 
24 06478 
47 03534 
44 06168 

22 07122 

23 07122 

32 00094 

33 03534 
26 07122 
21 00098 
26 07132 
21 07122 
26 04277 
16 064 74 
33 04986 
16 06462 
16 03953 
47 C5634 
26 07027 
73 06781 
26 0694 7 
39 06929 
16 03P65 
06 071 1 7 
15 06800 
44 04014 
33 071 1 5 
15 06800 
73 06819 

15 06822 

14 071 1 7 
47 04086 

16 071 1 7 
49 04122 
44 04122 
33 071 1 7 

15 06822 
73 06827 
44 04158 
49 05742 
44 04194 
26 06955 
49 05478 



OOORR 

01200 

OOORR 

01200 

07147 

00000 

0712K 

00000 

071 15 

07147 

000-1 

06476 

01200 

03534 

00893 

07157 

00000 

00000 

00098 

07217 

00098 

07137 

07122 

000-1 

00000 

0-000 

-7237 

00100 

06927 

06490 

06781 

00400 

-7237 

-7237 

OOOOJ 

071 15 

00000 

OOOOK 

071 15 

00001 

OOORR 

01200 

000-1 

00000 

071 17 

00000 

00002 

071 17 

05706 

00000 

05586 

06827 

OOOOO 



MOD 
OUT 



CM ADD6.99.6101 1 

BE A 14 

CM ADD 1 .99,6101 1 

BNE #+36 

TF ADD 1 . ADD6 

8 A14 

FSUB WORK +9. ADD 1 • 1 1 

SF A 15 

RNF *+24 , WORK+7 

TF ADD 1 « ADD6 

AM CT1,1,10 

C CT1.J1 

BNE A15 

BNF FIND-12.A15 

S ADD1.SADD+11 

M ADD1.S460 

SF 94 

CF A 1 5 

TF ADD 1,98 

A 98,T940 

TF ADD3.98 

A ADD 1 , ADD4 

TF PPADD+1 1 , ADD1 

TFM ROWCT.1,10 

CF DONE 

TFM CTCD,0,8 

TFM MOD+ 11, TAB 

BNC1 SW22 

TF PCD+98.BL+9R 

TNF TABL+18, TABCT 

TF PCD+18, TABL+18 

WACD PCD 

TFM 0UT+11.TAB.7 

TFL WORK+9. TAB, 7, OUTPUT CONVERSIONS 

TDM E14. 8-1, 1,11 

BNF *+36, WORK+7 

CF WORK+7 

TDM E14.8-l,2,ll 

TNF El 4. 8+ 18, WORK+7 

TDM E14, 8+21,1 

CM WORK+9, 99, 101 1 

BNE *+36 

TFM WORK+9, 01, 10 

B »+48 

BNF *+36,W0RK+9 

CF WORK+9 

TDM El 4. 8+2 1, 2 

TNF E14.8+26.W0RK+9 

BNF *+24,SW2 

B SW2+36 

BNF *+36, PRICE 

TF PCD+26,E14,8+26 

8 DON 
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04 194 


44 


04218 


05238 




BN C 


*+24 ,D0NF1 


04206 


49 


05418 


OOOCO 




B 


MAOUT 


04218 


44 


04254 


04986 




BNF 


*+36 tDONE 


04230 


26 


07027 


06827 




TF 


PCD+98.F14.3+26 


04242 


49 


05586 


OOOCO 




B 


PRICE 


04254 


26 


06955 


06827 




TF 


PCD+26.E14.8+26 


04266 


14 


03965 


-0000 


PPADD 


CM 


OUT+1 1 .0 


04278 


47 


0430? 


1 ?00 




BNF 


*+24 


042Q0 


26 


0698 1 


06667 




TF 


D CD+52 .PP+20 


04302 


1 1 


03965 


-00 10 




AM 


OUT+1 1.10 


04314 


39 


06929 


00400 




WACD 


PCD 


04326 


26 


070?7 


06927 




TF 


PCD+98.BL+9R 


"4338 


24 


03965 


06483 




C 


0UT+ 1 1 ,LA=T 


04350 


47 


03954 


1200 




BNE 


OUT 


04362 


1 1 


06474 


000-1 




AM 


ROWCT, 1.10 


04374 


73 


06645 


06474 




TNF 


R+l ,90WCT 


04386 


26 


06939 


06645 




TP 


PCD+1 n.P+1 


04398 


21 


06483 


07157 




A 


LAST . S460 


044 1 


21 


03953 


07157 




A 


MOD+1 1 ,S46C 


0442? 


1 1 


0646? 


OOO-l 




AM 


CTCD, 1.10 


04434 


24 


0646? 


06470 




C 


CTCn , J2 


04446 


47 


03930 


01200 




BNE 


MOD— 1 2 


04458 


44 


04986 


02658 




BNF 


DONE, XX 


04470 


12 


0713? 


-0010 


SWOUT 


SM 


ADD3, 10 


0448? 


06 


071 3K 


071?P 




TFL 


ADD3, ADD2.61 1 


04494 


22 


07127 


07157 




S 


ADD 2 , S460 


04506 


12 


0^132 


-0010 




SM 


ADD3.10 


04=1 8 


06 


071 3K 


0712P 




TFL 


ADD3 , ADD? .61 1 


04530 


1 1 


071 32 


-0020 




AM 


AD^ 3 ,20 


04542 


16 


06478 


000-0 




TFM 


CT1 ,0, 10 


04554 


1 1 


06464 


000-3 




AM 


CTPAR.3, 1 


04566 


?6 


064BB 


07 1 3? 




TF 


ADD , ADD3 


04578 


06 


071 1 7 


071 2K 




TFL 


WORK+9 . ADD 1 . 1 1 


04590 


09 


071 3K 


071 17 




FDI V 


ADD3,W0RK+9,6.ni VIDtS PIVOT 


04602 


1 1 


07132 


-00 10 




AM 


ADD3.10..ROW BY PIVOT POINT 


04614 


1 1 


06478 


000-1 




AM 


CT1 , 1 , 10 


046?6 


24 


06478 


06464 




r 


CT1 ,CTPAP 


04638 


47 


04590 


1 200 




BNE 


»-48 


04650 


26 


07152 


07217 




TF 


ADD7.T940. , VECTOR TRANSFORMATION 


0466? 


16 


06468 


000-0 




TFM 


CT2.0. 10 


04674 


26 


071 ?? 


07137 




TF 


ADD 1 . ADD4 


04686 


26 


07127 


07152 


MOD5 


TF 


ADD 2 . ADD 7 


04698 


16 


06478 


000-0 




TFM 


CT1 ,0. 10 


0471 


26 


071 ?.7 


071 ?? 




TF 


ADD4 . ADD 1 


04722 


?6 


071 3? 


06488 




TF 


ADD3 . ADD 


04734 


24 


06488 


07127 




C 


ADD . ADD2 


04746 


47 


04806 


01200 




BNE 


SW-24 


04758 


2 1 


071 ?? 


071 57 




A 


ADD 1 . S460 


04770 


21 


071 52 


07157 




A 


ADD7. C460 


04782 


21 


07127 


07157 




A 


ADD 2 . S460 


04794 


21 


07137 


07 157 




A 


ADD4 « S460 


048<">6 


06 


071 07 


071 3P 


/ 


T^L 


W0RK1+9, ADD4, 1 1 


0481 8 


06 


071 1 7 


071 "7 


TFL 


WORK + 9. WORK 1 +9 



04830 
04842 
04854 
04866 
04878 
04890 
04902 
04Q14 
04926 
04938 
04950 
ri4Q6? 
04974 
04986 
04998 
0501 
C5C22 
05034 
05C46 
05058 
05070 
05082 
05094 
051 06 
051 18 
05130 
05142 
05154 
05166 
05178 
05190 
0^202 
05214 
05226 
05238 
0^250 
05262 
05274 
05286 

05310 
05322 
05334 
05346 
05358 
05370 
05382 
05394 
05406 
054 18 
05430 
05442 
05454 



03 071 1 7 

02 0712P 

11 07127 

1 1 071 3? 

1 1 06478 

24 06478 

47 04818 

21 07152 

21 07122 

I 1 06468 
24 06468 
47 04686 
49 02178 
16 06478 
32 04986 
26 07147 
14 0714P 
47 05082 

II 07147 
43 061?4 
12 07147 
21 07147 
1 1 06478 
24 06478 
47 05022 
16 06478 
26 06488 
26 07122 
26 07027 
73 06781 
26 07003 
39 06"?9 
26 07027 
26 07027 
39 06929 
26 07027 
26 03965 
21 06488 
21 07177 
21 07167 

21 07172 
1 1 0647B 
24 06478 
47 03954 

22 06488 

32 05238 
26 07027 
26 03965 
49 03954 
26 07003 
26 06975 
39 06929 

33 0523e 



0713K 

071 17 

-0010 

-0010 

000-1 

06464 

1200 . 

07157 

071 57 

000-1 

06472 

01 200 

00000 

000-0 

00000 

O7?07 

000N0 

01200 

-001 1 

0714P 

-00 1 1 

07157 

000-1 

06476 

01200 

000-0 

07217 

07207 

06927 

06490 

06781 

00400 

06927 

06761 

00400 

06^27 

06488 

07157 

07157 

07! 57 

07157 

000-1 

06251 
01 200 
07157 
00000 
06927 
06488 
00000 
06827 
06799 
00400 
00000 
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SW FMUL WORK+9 . ADD3 .11 

FSUB ADDH . WORK+9 1 6 

AM ADD2.10 

AM ADD3.10 

AM CTl.1,10 

C CTl.CTPAR 

BNE SW-12 
A . ADD7.S460 

A ADD 1 .S460 

AM CT2.1.10 

C CT2.JM2 

BNE MODS 

B CJ 

DONE TFM CT1 .0.10, PUNCHES OUT 

SF DONE 

TF ADD6.T920 

NOSOL CM ADD6 . 50 » 6 1 

BNE *+48 

AM ADD6 ,11 

BD ERR . ADD6 .11 

SM ADD6.11 

A ADD6.S460 

AM CT1 , 1 , 10 

C CT1.J1 

BNE NOSOL 

TFM CTltO.10 

TF ADD.T940 

TF ADD1.T920 

TF PCD+98»BL+98 

TNF T ABL+ 18, TABCT 

TF PCD+74.TABL+18 

WACD PCD 

TF > PCD+98.BL+98 

TF PCD+98.BV+92 

DONE 1 WACD PCD 
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FINAL ANSWERS 



TF 
TF 



AM 
C 

BNE 

S 

SF 
TF 
TF 

B 

TF 
.TF 
WACD 
CF 



PCD+98.BL+98 
OUT+1 1 , ADD 
ADD.S460 
T470.S460 
T461 ,S460 
T462.S460 
CT1 , 1 , 10 
CT1 , J 
OUT 

ADD.S460 
DONE 1 

PCD+98.BL+98 
OUT+1 1 . ADD 
OUT 

PCD+74.E14.8+26 
PCD+46.MAX+14 

PCD 
DONE 1 
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05466 


49 


00402 


OOOOO 


B 


402 


05478 


1 5 


06467 


OOOOP 


DON TDM 


CT2— 1,7,11 


05490 


33 


05586 


OOOOO 


CF 


PR I CE 


05502 


1 4 


071 7P 


000-2 


D0N1 CM 


T470, 02,610 


055 1 4 


46 


05562 


1 200 


BE 


#+48 


05526 


25 


06468 


071 6P 


D0N2 TD 


CT2 , T46 1 ,11 


05538 


26 


06977 


06468 


TF 


PCD+48,CT2 


05550 


49 


05238 


OOOOO 


R 


DONE 1 


05562 


73 


06977 


071 IK 


DON3 TNF 


PCD+48.T462. 1 1 


05574 


49 


05238 


OOOOO 


B 


DONE 1 


05586 


26 


03965 


071 22 


PPICE TF 


OUT+ 1 1 , ADD 1 


05598 


2 1 


071 22 


071 57 


A 


ADD 1 , S460 


056 1 


32 


05586 


OOOOO 


SF 


PR I CE 


05622 


49 


03954 


OOOOO 


B 


out 


05634 


1 6 


06478 


000-0 


SW22 TFM 


CT 1 , , 1 , PUNCHES OUT ROW ZJ — CJ ONLY 


05646 


47 


04470 


00200 


BNC2 


SWOUT 


05658 


23 


0625 1 


071 57 


M 


J , S460 


05670 


32 


00095 


OOOOO 


SF 


95 


C5682 


2 1 


00099 


071 82 


A 


99 , T480 


05694 


26 


06488 


00099 


TF 


ADD , 99 


05706 


26 


03965 


06488 


SW2 TF 


OUT+ 11, ADD 


05718 


32 


05706 


OOOOO 


SF 


SW'2 


05730 


49 


03954 


OOOno 


3 


OUT 


05742 


73 


0678 1 


06490 


TNF 


TABL+18.TABCT 


05754 


34 


00000 


00 1 02 


RCT Y 




05766 


25 


06829 


00400 


TD 


E14«8+28,4Q0 


05778 


39 


06763 


00 100 


WATY 


TABL 


05790 


1 5 


06829 


OOOOO 


TDM 


E 1 4 • 8+28 , 


05802 


44 


04986 


02658 


BNF 


DONE , XX 


058 1 4 


33 


05706 


OOOOO 


CF 


SW2 


05R26 


49 


O4470 


OOOOO 


F? 


SWOUT 


05834 








DORG 


*— 3 


05863 




00030 




DC 


30,0 


05864 


32 


05836 


OOOOO 


PAP SF 


PAR— 28 


05876 


2' 


06464 


06464 


S 


CTPAP,CTPAR 


05888 


2 1 


06458 


071 62 


B A 


12, S470 , , SETS UP IDENTITY MATRIX 


05900 


06 


0645Q 


06441 


TFL 


I 2, ONE, 6 


05912 


1 1 


06478 


000-1 


AM 


CT1 ,1,10 


05924 


14 


05835 


OO0K4 


CM 


PAR-29,24, IC.DFTERMlKff WHETHFP 


05936 


46 


06044 


01200 


BE 


BACK, , .ARTIFICIAL VARIABLES 


05948 


14 


05835 


00OL3 


CM 


PAR-29,33, 10,ARt : . TO BE ADDED 


05960 


46 


06044 


01200 


BE 


BACK 


05972 


14 


05835 


000-4 


CM 


DAR-29,04, 10 


05984 


46 


06008 


01200 


BE 


*+24 


05996 


17 


06080 


-6497 


BTM 


ER, EREQ 


06008 


26 


06043 


06458 


TF 


*+35,I2 


06020 


12 


06043 


-0002 


SM 


*+23«2 


06032 


32 


0604L 


OOOOO 


SF 


*+l 1 , ,6 


06044 


26 


06279 


05863 


BACK TF 


IN+26,PAR-1 


06056 


32 


01698 


OOOOO 


SF 


SETT AB 


06068 


42 


OOOOO 


OOOOO 


q B 




06080 


34 


00000 


00102 


ER RCTY 




06092 


39 


0607R 


001 00 


WATY 


ER-1 , ,6 
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06104 


48 OOOOO 


OOOOO 




H 




061 16 


49 00402 


OOOOO 




B 


402 


06124 








DORG 


*-3 


061 24 


26 07027 


06927 


ERR 


TF 


PCD+98.BL+98 


061 36 


26 06967 


06633 




TF 


PCD+38 ,ERN0FE+38 


06 1 48 


39 06929 


00400 




WACD 


PCD 


06160 


49 00402 


OOOOO 




B 


402 


06168 








DORG' *-3 


06 1 68 


26 07142 


00893 






ADD5,SADD+1 1 


06180 


14 071 4K 


OOORR 


F I ND 


CM 


ADD5,99»6101 1 


061 92 


47 06228 


1 200 




BNE 


*+36 


06204 


11 07142 


-00 1 




AM 


Anns; - l n 


062 1 6 


49 06180 


OOOOO 








"6228 


26 0~M?2 


071 4? 






ADD 1 t ADD5 


06240 


49 03726 


OOOOO 




g F 


A 1 4+48 


06248 








DORG 


*-3 


06248 


00004 




I J 


DSS 




06253 


00080 






DAS 


80 


06249 


OOOOO 




I 




• I J+ 1 ♦ NUMBER OF VAR I ABLES XI* X2 ♦•••XI 


06251 


OOOOO 




j 


DS 


t i jTJtiNUj*ioc.K Ur KUWo 1 NC^ULJ I No MAXZ 


064 1 9 


00008 






DC 
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